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Applying  Particle  Tracking  Model 
in  the  Coastal  Modeling  System 

by  Honghai  Li,  Lihwa  Lin,  and  Mitchell  E.  Brown 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  describes  pro¬ 
cedures  to  set  up  a  particle  tracking  simulation  through  coupling  the  Coastal  Modeling  System 
(CMS)  and  Particle  Tracking  Model  (PTM)  in  the  Surface-water  Modeling  System  (SMS).  The 
step-by-step  CMS-PTM  setup  is  illustrated  by  a  coastal  application  at  Noyo,  CA. 

INTRODUCTION:  The  CMS  is  an  integrated  wave,  current,  sediment  transport,  and  morphol¬ 
ogy  change  model  for  coastal  and  inlet  applications  (Buttolph  et  al.  2006;  Lin  et  al.  2008)  that 
was  developed  by  the  Coastal  Inlets  Research  Program  (CIRP).  The  Particle  Tracking  Model 
(PTM)  was  developed  jointly  by  CIRP  and  the  Dredging  Operations  and  Engineering  Research 
(DOER)  program  as  a  tool  to  compute  the  fate  and  pathways  of  sediments  and  other  waterborne 
particulates  in  coastal  engineering  and  dredging  operations  (MacDonald  et  al.  2006;  Demirbilek 
et  al.  2008).  CMS-PTM  is  implemented  in  the  SMS  (Aquaveo  2010),  which  is  a  graphical  user 
interface  utility  for  PCs  developed  by  the  U.S.  Army  Corps  of  Engineers  (USACE). 

NOYO  NEARSHORE  SEDIMENT  PLACEMENT  PILOT  STUDY:  Noyo  Harbor  and  Noyo 
River  are  located  on  the  northern  California  coast.  Every  year  about  30,000-40,000  cy  of  beach- 
quality  sediment  is  dredged  from  the  Noyo  River  Federal  Navigation  Channel.  The  USACE  San 
Francisco  District  (SPN)  and  the  Coastal  and  Hydraulics  Laboratory  (CHL),  U.S.  Army  Engineer 
Research  and  Development  Center,  are  presently  collaborating  with  the  Noyo  Harbor  District 
and  stakeholders  in  a  pilot  study  to  investigate  a  potential  candidate  site  for  nearshore  placement 
of  dredged  sediment.  Upland  disposal  sites  used  in  the  past  will  soon  be  exceeded.  The  pilot 
study  proposes  a  dispersive  nearshore  site  3  to  8  miles  north  of  Noyo  Bay  along  the  northern 
California  coast  (Figure  1). 

To  evaluate  the  potential  impacts  of  nearshore  placement,  a  Short-Term  Fate  of  Dredged 
Material  Model  (STFATE)  and  CMS  were  applied  for  the  study.  SPN  has  applied  STFATE  to 
calculate  near-field  sediment  transport  and  suspended  sediment  concentration  during  the  place¬ 
ment  of  the  dredged  material.  The  USACE  CIRP  applied  the  CMS  and  CMS-PTM  to  simulate 
sediment  transport  for  combined  waves  and  currents  under  various  ambient  conditions  and 
project  alternatives.  The  model  results  will  provide  the  technical  information  necessary  to  eval¬ 
uate  a  demonstration  site  that  is  economically  feasible  for  the  optimum  sediment  placement 
location. 

CMS-PTM:  The  CMS  calculates  water  levels,  currents  and  waves  through  the  coupling  of  a 
hydrodynamic  model,  CMS-Flow  and  a  wave  spectral  model,  CMS-Wave.  CMS-Flow  also 
simulates  sediment  and  salinity  transport,  and  morphology  change. 
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Figure  1.  CMS  domain,  grid,  and  bathymetry. 


CMS-Flow  is  driven  by  time-dependent  water  surface  elevation  at  the  offshore  open  boundaries, 
and  wind  forcing  over  the  surface  boundary.  Directional  wave  data  for  input  to  CMS-Wave  are 
available  from  the  National  Data  Buoy  Center  (NDBC,  http://www.ndbc.noaa.gov)  and  Coastal 
Data  Information  Program  (CDIP,  http://cdip.ucsd.edu).  The  wave  conditions  are  specified  at  the 
seaward  boundaries. 
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CMS-Flow  requires  three  input  files  and  CMS- Wave  requires  six  input  files  through  the  SMS. 
Each  of  the  CMS-Flow  and  CMS-Wave  simulations  will  generate  its  own  output  file  (Figure  2). 


Figure  2.  CMS  files. 


The  CMS-Flow  input  and  output  files  are  required,  while  the  CMS-Wave  output  files  are 
optional  for  setting  up  a  CMS-PTM  simulation.  After  loading  the  required  hydrodynamic  and/or 
wave  files,  users  can  follow  the  steps  below  to  set  up  the  CMS-PTM  simulation  for  the  Noyo 
application. 

1 .  Native  Sediment  Dataset:  CMS-PTM  requires  spatially  varying  native  bed  sediment  data 
specified  as  datasets,  including  the  grain  size  information  of  d35,  d50,  and  d90  for  each  cell, 
in  the  CMS-Flow  grid. 

i)  Select  and  open  the  CMS-Flow  control  file,  “flow-exist-Alt2.cmcards  ”  (Figure  2). 

ii)  Choose  Data  \  Data  Calculator  in  the  CMS-Flow  menu. 

iii)  Name  the  first  dataset,  d35,  in  the  Output  data  set  name  and  specify  a  constant  value, 
0. 1  mm  in  the  Calculator. 

iv)  Click  Compute  to  generate  the  dataset. 

v)  Repeat  steps  iii)  and  iv)  to  generate  the  datasets  for  d50  and  d90  (Figure  3),  (for  the 
Noyo  simulation,  d50=0.2  mm  and  d90=0.4  mm). 
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Figure  3.  Generate  the  datasets,  d35,  d50,  and  d90. 

vi)  Click  Done  to  exit  the  Data  Calculator. 

vii)  Highlight  the  datasets,  d35,  d50,  and  d90,  right-click,  and  select  Export  Datasets  (Fig¬ 
ure  4). 

viii)  Select  XMDF  file  in  the  File  Type,  check  All  time-steps  in  the  Time-steps,  and  click  the 
Filename  to  save  the  file  “NativeSediments.h5”  in  the  present  working  folder. 

2.  Land-water  Interface  Boundary  Strings:  CMS-PTM  requires  specifications  of  both 
open  and  closed  boundaries.  Open  boundaries  are  specified  in  CMS-Flow,  so  closed  bounda¬ 
ries  need  to  be  added  for  CMS-PTM  (future  versions  of  SMS  will  automatically  generate  the 
closed  boundaries).  The  steps  to  perform  this  task  are  described  as  follows: 

i)  Choose  Cellstring  |  Generate  Along  Boundary  (this  step  will  generate  redundant  boun¬ 
dary  strings  overlapped  with  existing  strings  along  open  boundaries)  (Figure  5). 

ii)  Choose  the  Select  Cellstring  S  tool  and  select  the  Land  boundary  along  the  open 
boundary. 

iii)  Right-click  and  choose  Delete  (Figure  5). 
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Figure  4.  Generate  the  native  sediment  file. 


iv)  Repeat  steps  ii)  and  iii)  to  delete  other  Land  boundary  cellstrings  overlapped  with 
existing  open  boundaries  if  more  than  one  open  boundaries  exist  in  CMS-Flow. 

v)  Save  the  project. 


3. 


Sediment  Sources:  Particles  can  be  released  from  different  source  options.  In  CMS-PTM, 
sediment  sources  can  be  specified  as  points,  lines,  or  areas.  For  the  Noyo  application,  an  area 
source  was  designated  to  represent  the  nearshore  placement  of  dredged  materials  from  a 
barge.  The  steps  for  creating  the  source  file  in  SMS  are: 


i) 


Switch  to  the  Map  Module  by  clicking  the  icon 


ii)  Right-click  the  default  coverage  and  choose  Type  |  Models  \  PTM  after  generating  the 
rectangular  area  for  the  sediment  release  (Figure  6). 

iii)  Choose  the  Select  Feature  Polygon  icon  ♦***  and  click  on  the  area  source. 


iv)  Right-click  the  polygon,  select  Attributes ,  and  the  Feature  Object  Attributes  will  pop 
up. 


v)  Select  the  Horizontal  Polygon  Source  in  the  Type  (Figure  7). 

vi)  Specify  Date/Time  for  the  particle  release,  Parcel  Mass  (kg),  released  mass  Rate 
(kg/m2/s),  and  Median  Grain  Size  (mm)  of  released  particles  (Figure  7). 

vii)  Save  the  area  source  specification  as  a  map  file. 
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Figure  6.  Link  source  file  to  CMS-PTM. 


Figure  7.  Define  particle  release  scenario  from  the  area  source. 


4. 


CMS-PTM  Setup:  After  the  hydrodynamic  data,  native  sediment  data,  source  file,  and 
boundary  cellstrings  are  generated,  users  can  follow  the  steps  below  to  specify  the  CMS- 
PTM  parameters  and  set  up  the  CMS-PTM  simulation. 


i) 


Switch  to  Particle  Module  by  clicking  the  icon 
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ii)  Choose  PTM  \  New  Simulation. 

iii)  Choose  PTM  \  Model  Control. 

•  Simulation  time  (Figure  8). 

i)  Specify  the  starting  time  of  the  completed  CMS-Flow  simulation  in  the  Hydrodynam¬ 
ics. 

ii)  Specify  the  starting  and  ending  times  of  a  CMS-PTM  run  in  the  Simulation. 

iii)  Specify  the  time-step  for  the  CMS-PTM  simulation  and  the  output  time  interval  in  the 
Particle  Output. 

iv)  Specify  how  often  to  update  the  sediment  parameters  and  the  hydrodynamic  variables. 

The  starting  time  of  CMS-Wave  and  the  time-step  for  the  wave  forcing  update  need  to  be  speci¬ 
fied  if  waves  are  included  in  the  simulation  and  selected  in  the  Files  page. 


Figure  8.  CMS-PTM  model  control  pages  (Time). 

•  Links  to  hydrodynamic  and  sediment  data  (Figure  9). 

i)  Click  Select  paths  under  the  XMDF  Path  and  set  Alt2/Datasets/Depth  for  the  Depth 
Path. 

ii)  Select  the  CMS  solution  file,  flow-exist-alt2_sol.h5,  in  the  working  folder  for  the 
Hydrodynamics. 
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Figure  9.  CMS-PTM  model  control  pages  (Files). 

iii)  Set  Simulation/Current _Velocity  for  the  Velocity,  and  Simulation/  Water _Elevation  for 
the  Water  surface  elevation. 

iv)  Select  default  coverage  in  the  Sediment  source  and  name  the  part. source  in  Select  file¬ 
name. 

v)  Select  create  when  model  is  run  in  the  Neighbor  and  name  the  part. neighbors  in  Select 
filename. 

vi)  Include  Waves  and  breaking  by  selecting  CMS- Wave  and  click  options  to  obtain  wave 
input  for  CMS-PTM  (Figure  10). 

vii)  Click  the  icon  ^1  to  browse  the  working  folder  and  to  link  the  wave  simulation  file, 
wave-exist-Alt2  .sim,  under  the  Parent  Grid  Geometry  \  CMS  simulation  file. 

viii)  Click  the  icon  Load  ^ ,wav2  I  to  link  the  wave  input  file,  wave-exist-Alt2.wav,  under  the 
Parent  Grid  Solution  Data. 

ix)  Specify  the  Number  of  Spectra  in  file. 

x)  Click  OK  to  complete  the  specifications  of  the  page. 

Because  the  wave  simulation  is  included  in  CMS-PTM,  it  is  necessary  to  return  to  the  previous 
(Time)  page  to  specify  the  Start  time  and  the  Time-step  under  the  Waves  (Figure  8). 
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Figure  10.  Specifications  of  wave  input. 

•  Computation  parameters  and  output 

Most  of  the  options  in  the  Computations  page  are  default  for  this  Noyo  application,  except  that 
the  3D  mode  is  selected  (Figure  11).  In  the  Output  page  (Figure  12),  Print  default  keywords  is 
selected  in  the  File  Output,  and  Source  and  State  are  selected  in  the  Particle  Output. 

5.  CMS-PTM  Results:  As  described  in  the  previous  sections  for  the  Noyo  application,  CMS- 
PTM  was  designed  to  investigate  the  fate  and  transport  of  dredged  materials  under  the  com¬ 
bined  influence  of  waves  and  currents  from  proposed  placement  sites. 

The  placement  site,  a  rectangular  area  of  86,600  yd  ,  is  located  about  4  miles  north  of  Noyo  Har¬ 
bor  near  MacKerricher  State  Park.  The  CMS-PTM  simulation  is  for  August  1-30,  2008.  In  this 
application,  clay  materials  were  released  at  a  level  of  1  m  above  the  ocean  bottom  at  a  rate  of 
twice  daily  for  12  days.  Each  release  lasted  10  minutes  to  mimic  the  barge  placement  operation. 
Figures  13  and  14  show  the  snapshot  of  calculated  clay  particle  distributions  at  the  end  of  the 
release  (day  12)  and  the  end  of  simulation  (day  30),  respectively.  The  simulation  shows  that  clay 
particles  either  follow  the  waves  propagating  towards  shore  or  move  southward,  driven  by  the 
northerly  wind  in  this  summer  period.  A  good  portion  of  the  fine  sediments  can  move  past  Noyo 
Bay  to  the  model’s  south  boundary  and  some  out  of  the  model  domain  through  the  simulation.  At 
the  end  of  the  simulation,  about  65  percent  of  the  released  clay  particles  are  considered  “dead,” 
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which  means  that  they  are  either  permanently  buried  at  the  sea  bed  or  have  moved  out  of  the 
model  domain. 


Figure  1 1 .  Specifications  of  CMS-PTM  parameters. 


Figure  12.  Selections  of  CMS-PTM  output. 
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Figure  13.  An  oblique  view  of  calculated  clay  particle  distribution  at  the  end  of  the  release  (day  12). 
The  white  square  indicates  the  release  area. 


Figure  14.  An  oblique  view  of  calculated  clay  particle  distribution  at  the  end  of  the  simulation  (day  30). 
The  white  square  indicates  the  release  area. 
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ADDITIONAL  INFORMATION:  This  CHETN  was  prepared  and  funded  under  the  Coastal 
Inlets  Research  Program  (CIRP)  and  was  written  by  Dr.  Honghai  Li  ( Honghai.Li@usace.army . 
mil,  voice:  601-634-2840,  fax:  601-634-3080),  Dr.  Lihwa  Lin  ( Lihwa.Lin@usace.army.mil ),  and 
Mitchell  E.  Brown  ( Mitchell.E.Brown@usace.army.mil )  of  the  U.S.  Army  Engineer  Research 
and  Development  Center  (ERDC),  Coastal  and  Hydraulics  Laboratory  (CHL).  The  CIRP  Pro¬ 
gram  Manager,  Dr.  Julie  D.  Rosati  ( Julie.D.Rosati@usace.army.mil ),  the  assistant  Program  Man¬ 
ager,  Dr.  Nicholas  C.  Kraus,  and  the  Chief  of  the  Coastal  Engineering  Branch  at  CHL, 
Dr.  Jeffrey  P.  Waters  reviewed  this  CHETN.  Files  for  the  example  may  be  obtained  by  contact¬ 
ing  Dr.  Li.  This  CHETN  should  be  cited  as  follows: 

Li,  H.,  L.  Lin,  and  M.  E.  Brown.  2011.  Applications  of  Particle  Tracking  Model 
in  the  Coastal  Modeling  System.  Coastal  and  Hydraulics  Engineering  Technical 
Note  ERDC/CHL  CHETN-IV-78.  Vicksburg,  MS:  U.S.  Army  Engineer  Research 
and  Development  Center.  An  electronic  copy  of  this  CHETN  is  available  from 

http: //chi.  erdc.  usace. army.mil/chetn  or  http.V/cirp.  usace.  army. mil/pubs/chetns. html. 
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